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ABSTRACT 

In a previous paper (Oiling 1995, AJ, 110, 591) a method was developed to determine 
the shapes of dark matter halos of spiral galaxies from an accurate determination of 
the rotation curve, the flaring of the gas layer and the velocity dispersion in the H I . 
Here this method is applied to the almost edge-on Scd galaxy NGC 4244 for which the 
necessary parameters are determined in the accompanying paper (A J, Aug. 1996). 

The observed flaring of the H I beyond the optical disk puts significant constraints on 
the shape of the dark matter halo, which are almost independent of the stellar mass-to- 
light ratio. NGC 4244's dark matter halo is found to be highly flattened with a shortest- 
to-longest axis ratio of 0.21q ?■ If the dark matter is disk-like, the data presented in this 
paper imply that the vertical velocity dispersion of the dark matter must be 10% - 30% 
larger than the measured tangential dispersion in the H I . Alternatively, the measured 
flaring curve is consistent with a round halo if the gaseous velocity dispersion ellipsoid 
is anisotropic. In that case the vertical dispersion of the gas is 50 - 70% of the measured 
tangential velocity dispersion. 
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1. INTRODUCTION 

Although rotation curves of spiral galaxies have 
been used as evidence for the presence of dark mat- 
ter (DM), little is known about the nature, extent 
and actual distribution of the DM halos (e.g., van 
Albada et al. 1985; Lake & Feinswog 1989). As mea- 
surements of the equatorial rotation curve probe the 
potential in only one direction, they provide no infor- 
mation about the shape of the DM halos. 

In a previous paper, Paper I (Oiling 1995), a 
method was developed to determine the shape of the 
dark matter halo from the gaseous velocity disper- 
sion and the radial variation of the thickness of gas 
layers (flaring). This is accomplished by comparing 
the measured flaring with that expected from a self- 
gravitating gaseous disk in the axisymmetric poten- 
tial due to the stellar disk and (flattened) DM halo. 
The method of calculation which produces the most 
reliable model flaring curve is the global approach, in 
which the potential is calculated from the total mass 
distribution of the galaxy (Paper I) . In this paper the 
method is applied to the galaxy NGC 4244 for which 
the basic parameters were determined in the accom- 
panying paper (Oiling 1996, hereafter Paper II). 

In the past, several methods have been used to 
determine the shapes of dark matter halos. Hofncr 
& Sparke (1994) analyzed the warping behavior of 
H I disks and concluded that only one (NGC 2903) 
of the five systems studied requires a DM halo as 
flattened as E4Q On the other hand, in studies of 
polar ring galaxies (Sackett & Sparke 1990; Sackett 
et al. 1994; Sackett & Pogge 1995) substantially flat- 
tened DM halos are found (E6-E7 for NGC 4650A, 
E5 for A0136-0801). The shape of the dark halo of 
the Milky Way has been estimated (E0 - E7) from the 
kinematics of extreme Population II stars (Binney & 
Ostriker 1987b; Sommer-Larsen & Zhen 1990; van der 
Marel 1991; Amendt & Cuddeford 1994). Steiman- 
Cameron et al. (1992) used the dynamics of the pro- 
cessing dusty disk of the SO galaxy NGC 4753 to infer 
the shape of its DM halo (El). 



Van der Kruit (1981) pioneered the use of flaring 
measurements to determine the mass of stellar disks, 
and found that the scale length of the total matter 
(luminous plus dark) was similar to the scale length 
of the light distribution and concluded that the mass- 
to-light ratio does not vary significantly with radius. 
Since mass models where all the mass is concentrated 
in the disk require large radial changes in the mass- 
to-light ratio {M/C ) to reproduce the observed flat 
rotation curves, the DM halo is not as highly flat- 
tened as the stellar disk. A similar conclusion was 
reached for NGC 4565, NGC 891 (Rupen 1991), and 
the Milky Way (Knapp 1987; Merrifield 1992; Mal- 
hotra 1994, 1995). As we can measure the light dis- 
tribution of the stellar disk only and have no prior 
knowledge of the mass-to-light ratio of stellar disks, 
the relative contributions of luminous and dark mat- 
ter to a given rotation curve are not known. Models 
of galaxy formation which include simple estimates 
of the star formation history (e.g., Larson & Tinsley 
1978) predict M/C 's which are in the range of val- 
ues derived from observed rotation curves (e.g., Rubin 
et al. 1985; Athanassoula et al. 1987; Bottema 1988; 
Broeils 1992, Chap. 12), but different star formation 
histories can easily result in mass-to-light ratios which 
differ by a factor of a few (e.g., Tinsley 1981; Worthey 
1994; de Jong 1995, Chap. 4). An upper limit to 
the mass-to-light ratio, and hence a lower limit to the 
amount of DM, can be obtained by assuming that the 
peak of the observed rotation curve is due to the stel- 
lar disk only, the so called maximum-disk hypothesis 
(van Albada & Sancisi 1986). Other M/C estima- 
tors have been used (van der Kruit 1981; Efstathiou 
et al. 1982; Athanassoula et al. 1987; Martinet 1988; 
Bottema 1993), typically yielding a stellar disk 50% - 
100% as massive as the "maximum-disk". Employing 
the maximum-disk model, Broeils (1992, 1995) finds 
that the ratio of dark to luminous matter lies between 
0.5 and 10, where the larger values are found in low 
mass systems. Unfortunately, the flaring data pre- 
sented in this paper cannot be used to constrain the 



1 A spheroidal system with shortest-to-longest axis ratio 
c/a (q p ) of the density contours has a shape En, with 
n such that q p — 1 — n/10 



mass of the stellar disk of NGC 4244 (Sec. ^3j) be- 
yond the limits imposed by the shape of the rotation 
curve. 

In a study similar to the present one, Rupen (1991) 
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analyzed high resolution H I observations of two edge- 
on galaxies, NGC 4565 and NGC 891. For these 
galaxies he finds that the flaring data imply stellar 
disks with masses two-third and one-half as massive 
as a "maximum-disk" , consistent with van der Kruit's 
(1981, 1988) results for NGC 891. In addition, Ru- 
pen found that non-thermal pressure gradients 
need probably not be included in the equation of hy- 
drostatic equilibrium ([[]). The fact that NGC 891 
and NGC 4565 have similar luminous mass distri- 
butions, similar gas layer widths and similar ther- 
mal pressure gradients (^ k ) but are estimated to 
have very different non-thermal pressure gradients 

( rffiit / rf-Pth I r~ I ■ dP nt I rffth I r~ A\ ar _ 

v dz I dz INGC4565 ~ 2 ' dz I dz INGC891 ~ ' 

gues against the importance of non-thermal pressure 
gradients in balancing the vertical force. 

Below I summarize the relevant properties of NGC 
4244 as determined in Paper II and references therein. 
The current starformation rate in NGC 4244 is low. 
Apart from a region between 6.5 and 8 kpc to the 
south-west, the radial light profile is well represented 
by an exponential distribution: extinction is probably 
not very important. The properties of the neural hy- 
drogen layer were derived from sensitive high resolu- 
tion VLA0 observations. The total H I distribution 
is fairly regular and shows a gentle warp. Several 
strong H I surface density variations are seen in the 
inner disk. Beyond the optical disk the surface den- 
sity on both sides of the galaxy is similar, validating 
the use of the average surface density profile in the 
model calculations. The same holds for the gaseous 
velocity dispersion, and the inferred inclination an- 
gle. Superimposed upon a gradual increase of the 
gas layer width, the derived thickness of the gas layer 
shows a characteristic increase/decrease on the in- 
side/outside of spiral arms. These features are prob- 
ably not real but are likely an artifact (induced by 
spiral arm streaming motions) of the thickness-and- 
inclination derivation technique (Paper II) . NGC 4244 
is an intermediate mass system with a rotation speed 
of about 100 km s _1 . The large inclination (~ 84° 5) 

2 The VLA of the National Radio Astronomy Observatory 
is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, 
Inc. 



of NGC 4244 in combination with the symmetric sur- 
face density distribution, allowed for the accurate de- 
termination of the rotation curve. NGC 4244 is one 
of the few galaxies to show, on both sides, the char- 
acteristic drop in rotation velocity beyond the optical 
disk. This decline in rotation speed is more typical 
of luminous fast rotators (Casertano & van Gorkom 
1991). The slow rise in the inner part and the decline 
of the rotation curve beyond the optical disk limit the 
mass-to-light ratio of NGC 4244 to 50% - 100% of the 
maximum-disk value. Due to the limited extend of the 
measured rotation curve the dark halo's core radius 
(R c ) and central density (po) could be determined to 
within a factor of 50 and 900 only, respectively. This 
large uncertainty in the dark-halo's density distribu- 
tion does not preclude the determination of the halo's 
shape (Paper I). 

In Sec. |]] I describe how, in conjunction with self- 
consistent mass models, the observed flaring can be 
used to constrain the shape of the dark matter halo. 
The shape of the halo is then determined in Sec. |2.1.| . 
Some caveats related to various step along the way are 
discussed in Sec. ED. The implications of the highly 
flattened halo for the amount and nature of dark mat- 
ter are discussed in Sec. [h| 

2. THE MASS MODEL 

The dependence of the thickness of the gas layer 
upon the shape of the dark matter halo is described 
in Paper I. Here I will only present four of the most 
salient features. That the thickness of a gas layer de- 
pends on the axial ratio (q p = c/a) of a dark halo 
can be easily understood. For an axisymmetric ellip- 
soidal dark matter distribution with fixed R c and po 
[likeEq. (§], the mass enclosed by an ellipsoid of ax- 
ial ratio q p , and hence the rotation speed, decreases 
with decreasing q p . Thus, for a given observed ro- 
tation curve, the densities of flattened halos have to 
be larger than for a round halo . The vertical force, 
proportional to the mass enclosed (Gauss' theorem), 
increases as well so that flattened halos have thinner 
gas layers. The second important result obtained in 
Paper I is that the mass-to-light ratio (or equivalently, 
the mass) of the stellar disk has virtually no influence 
on the thickness of the gas layer beyond the stellar 
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disk. Given the fact that the dark matter density 
distribution is poorly known (van Albada et al. 1985; 
Lake & Feinswog 1989; Paper I; Paper II) this may 
seem surprising at first. However, close to the plane 
the vertical force is roughly proportional to the ra- 
dial forced which is, by construction, the same for all 
disk-halo combinations that reproduce the observed 
rotation curve: beyond the optical disk the depen- 
dence on halo flattening dominates. Thirdly, the self- 
gravity of the gas can contribute significantly to the 
vertical force beyond 2-3 optical scale lengths. And 
finally, the model gas layer widths are best calculated 
using the global approach, in which the vertical force 
is determined from the global mass distribution of the 
galaxy. 

Comparing the thickness measurements beyond 
the optical disk with model flaring curves, calculated 
for a series of mass models with varying halo flatten- 
ing, yields the flattening of the DM halo. The basic 
assumptions are that the vertical force K z (z) is bal- 
anced by gradients in the thermal pressure (Pth) and 
that the gaseous velocity dispersion ellipsoid is round 
and independent of height above the plane, so that 
the vertical Jeans equation [e.g., Paper I, Eq. (1)] 
can be approximated by the equation of hydrostatic 
equilibrium: 



et al. 1994, respectively). Here a constant vertical 
scale height is assumed. For the vertical distribution 
I use an exponential as well as a sech 2 (z/zo) distri- 
bution. The effects of a different choice for the ver- 
tical stellar distribution are discussed in Sec. 3.3.1 



To avoid extinction features (Paper II), a fit to the 
stellar radial surface brightness distribution is used: 
L{R) = L(0)exp(-R/h R ) for R < 10 kpc, followed 
by a linear decline till 11 kpc and L(R) — be- 
yond 11 kpc. A constant mass-to-light ratio (M/Cb ) 
is assumed. The stellar density distribution is thus 
given by the product of the mass-to-light ratio and 
the adopted radial and vertical distributions. 

The gaseous surface density used is the sum of the 
H I surface density distribution, a crude estimate of 
the H2 surface density distribution, and 25% He by 
mass (Paper II). The volume density is found from 
£ gas (i?), (r(-R)gas and the total potential, via the equa- 
tion of hydrostatic equilibrium. Details of this proce- 
dure can be found in Paper I. A Gaussian fit to this 
density distribution yields the model gas layer widths 
(but note that deviations from the Gaussian distri- 
bution occur inside the optical disk, see footnote 4 



th 



dz 



Pgas -^z 



(1) 



in Paper I). In Sec. 2.1. these model widths will be 
compared to the measured flaring curve. 

The dark halo volume density is parameterized as 
a flattened isothermal sphere with core radius R c {q p ), 
central density po(q p ), and flattening q p (= c/a) (e.g., 
Sackett & Sparke 1990; Sackett et al. 1994; Paper I): 



Input to these mass models are the gaseous velocity 
dispersion, the observed rotation curve, and the den- 
sity distributions of the dark and luminous matter. 
Below I describe how to arrive at these density distri- 
butions. 

Optical studies of edge-on spiral galaxies indi- 
cate that the vertical scale height of the stellar disk 
is constant as a function of galactocentric distance 
(van dcr Kruit & Searle 1981a,b 1982a,b; Shaw & 
Gilmore 1990; Barteldrees & Dettmar 1994), although 
there are theoretical and observational indications 
that the shape of the vertical distribution may change 
with radius (Burkert et al. 1992 and van Dokkum 

3 K z =z/Vz 2 +R 2 F tot « z/R oc -§ V D 2 bs 



p(R,z;q p ) 



Po(q P ) Rcilp 



R 2 M + R 2 + {z/q p f 



(2) 



The dependency of R c and po upon q p is such that 
the rotation curve of any flattened density distribu- 
tion described by Eq. (|J) is practically indistinguish- 
able from its round equivalent (Paper I). In the flat 
rotation curve regime, the thickness of the gas layer 
in a density distribution like Eq. (0) is given by: 



FWHM Z (R) « 2.35 



2.4 q 
lA + q 



R 2 + R, 



2 

c J 



(3) 
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with Vrot.max the asymptotic rotation speed of the DM 
halo [Paper I, Eq. (D31)]. 

As the mass models described above are fully spec- 
ified by two independent free parameters, the stellar 
mass-to-light ratio and the flattening of the DM halo, 
the flaring behavior can be investigated as a func- 
tion of these parameters. For every choice of the 
stellar M/C , the DM core radius and central den- 
sity are uniquely determined by two points on the 
observed rotation curve (e.g., at 2.3/ir and 8/ir, see 
Paper I, Appendix B). However, the predicted flar- 
ing beyond the optical disk is essentially independent 
of the choice of M/C (Paper I; Fig. 0). The shape 
of the DM halo is not a parameter in such rotation 
curve decompositions, although different halo Batten- 
ings require somewhat different values of R c and q p 
[Paper I, Eqs. (B8) and (B9)]. In Paper II it was 
found that the DM-halo parameters are not well con- 
strained: possible core radii range from 0.3 to 12 kpc 
and central densities between 1.3 and 0.001 M Q /pc 3 
are allowed; small M/C 's correspond to small core 
radii and large halo densities. Again, notwithstand- 
ing this enormous allowed range in halo parameters, 
the thickness of the gas layer in such disk-halo models 
(with the same flattening) is very similar. 

The Dark Matter halo model [Eqn. (§)] is not 
unique (e.g., van Albada et al. 1985; Lake & Feinswog 
1989; Navarro, Frenk & White 1996). However, all 
DM density distributions with similar rotation curves 
share the feature that flatter distributions have larger 
vertical forces and hence thinner gas layers. Thus, the 
flattening determined in this paper is an indicator of 
the true flattening of the true DM density distribu- 
tion. In the discussion (Sec. I investigate the ex- 
treme possibility that the dark matter is disk-like. 

2.1. The Shape of the Dark Halo 

Now we can compare the calculated gas scale 
heights for several mass models with various M/C 
and q p with the observed widths (Fig. |l|). The trun- 
cation of the stellar disk, the gaseous self-gravity and 
the measured gaseous velocity dispersions have been 
included in the model calculations. Inspecting this 
figure we notice general correspondence between the 
observations and the model flaring curves, except for 



the 7 kpc hump, where the thickness of the gas layer 
increases by a factor 2.7 within 1 kpc (on both sides 
of the galaxy) . Similar features are seen in the Milky 
Way (Knapp 1987) and NGC 891 (Rupen 1991) but 
are absent in NGC 4565 (Rupen 1991). The steep 
thickness gradients in NGC 891 and NGC 4244 seem 
to occur just inside and just outside H I surface den- 
sity enhancements. Such sudden changes in apparent 
thickness can arise from streaming motions induced 
by spiral structure (Paper II, Appendix C). 

For the best-fit halo flattening, the maximum-disk 
and the minimum-disk model curves seem to repre- 
sent the flaring data equally well (upper and lower 
panel, respectively). This is corroborated by the fact 
that the model with the largest reduced \ 2 li es within 
the 1-ct bound (i.e., \ 2 < Xmin + -0 °f the model with 
the smallest reduced Xmin- The north-eastern, south- 
western, and the average flaring curves all yield the 
same estimate for the mass-to-light ratio. Thus, for 
NGC 4244, the measured flaring does not constrain 
the mass-to-light ratio of the stellar disk beyond the 
limits set by the shape of the rotation curve. 

On the other hand, the flattening of the halo is 

much better constrained. For all stellar disk masses, 

the minimum x 2 is found at the same halo flattening, 

with the same formal uncertainty. The uncertainty 

in the inclination determination in the outer parts of 

the galaxy (Paper II) has only a small effect on the 

inferred DM-halo shape: for the varying-inclination 

case I find g p =0.2 ± 0.1, the fixed-inclination case 

yields g p =0.3j^i- There are several reasons why I 

believe that the small decrease in inclination beyond 

the optical disk is probably real: 1) it is seen on both 

sides of the galaxy, 2) it is seen (with lower S/N) when 

using the low velocity channels (8-14) only, as well as 

in the "edge" channels (15-20) only. Nevertheless, I 

include the fixed-inclination case in the final solution: 

a - 2 +0 - 3 
q p — u.^_ v 

3. CAVEATS 

We have seen that the dark halo of NGC 4244 is 
significantly flattened. Below I discuss the main un- 
certainties in this result. 
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3.1. The Gaseous Velocity Dispersion 

As the calculated thickness of the gas layer and the 
inferred halo flattening depend, respectively, linearly 
and inversely quadratically on the velocity dispersion 
of the gas^, the velocity dispersion of the gas is an 
important parameter in the analysis presented above. 
This strong dependence in illustrated in Fig. ^, where 
I compare model gas layer widths calculated for the 
lower and upper bounds to the gaseous velocity dis- 
persion found in other, more face-on, systems (van der 
Kruit & Shostak 1982; Dickey et al. 1990; see Kam- 
phuis 1993, Chap. 12, for a review). It is striking 
that the shape of the halo cannot be determined if no 
velocity dispersion measurements are available: mea- 
suring the gaseous velocity dispersion is instrumental 
in the determination of the shape of the dark matter 
halo. From Fig. ^ we can see that forcing the flaring 
measurements to be consistent with a round DM-halo 
model requires a radial decline of the velocity disper- 
sion which is inconsistent with the dispersion mea- 
surements. Given that line-of-sight smearing effects 
are unimportant and that the data from which the 
tangential velocity dispersion was determined has a 
high signal to noise ratio (S/N > 10 for R < 14 kpc; 
Paper II), it is highly unlikely that the measured ve- 
locity dispersion curve is consistent with the velocity 
dispersion needed to infer a round halo. 

The flaring data can be consistent with a round 
halo if the vertical dispersion (er zz ) does not equal 
the measured tangential dispersion (age). The mea- 
sured tangential dispersion curve can thus be con- 
sistent with a round DM halo if (<J zz /crgg) decreases 
radially (down to 0.55 at the last measured point). 
However, it is unclear whether such a large anisotropy 
in the velocity dispersion of the H I beyond the opti- 
cal disk is consistent with the observations. For the 
molecular component in the solar neighborhood this 
ratio is not well established: cr zz = 5.7 ±1.2 for high 
latitude CO clouds (Blitz et al. 1984) while determi- 
nations of agg from tangent point measurements vary 
between 3 ± 0.7 (Clemens 1985) and 7.8 ± 4 (Malho- 
tra 1994). 

Assuming that an anisotropy in the gaseous ve- 
4 From Eq. (||) it follows that FWHMl oc o^ as q p . 



locity dispersion will be reflected in the motions of 
young stars formed from this gas, proper motion stud- 
ies of OB associations may be informative. On very 
small scales (~ 0.5 pc) there is some evidence that 
the velocity dispersion ellipsoid in the Orion Nebula 
Cluster is anisotropic ((Tft/cr; = 2.2/2.5 = 0.9 ± 0.1), 
while it is uncertain how much this value is affected 
by the expansion of the cluster (Jones & Walker 
1988). On larger scales (~ 20 pc), the anisotropy 
of a proper motion selected sample of stars is simi- 
lar (2.9/3.2=0.9 ± 0.3; McNamara et al. 1989), while 
a large proper motion and distance selected sample 
(Warren & Hesser 1977) covering a large volume (175 3 
pc 3 ) centered on Orion OBI shows an increase of dis- 
persion values but no evidence for a velocity disper- 
sion anisotropy (17.3/16.7 = 1.04 ± 0.2). Mihalas & 
Binney (1981, p. 423) report a small anisotropy in 
the dispersion for 0-B5 stars (agg/cr zz = 11/9 = 1.2) 
with an unknown error. In summary, the velocity 
dispersion ellipsoid of young OB stars, and presum- 
ably of the interstellar medium, is not significantly 
anisotropic in the solar neighborhood. 

Possibly the best way to determine the velocity 
dispersion ellipsoid, is by comparing the velocity dis- 
persions in a sample of face-on and edge-on galaxies. 
We (van Gorkom, Rupen & Oiling) are in the process 
of doing so. 

3.2. The Determination of the Inclina- 
tion 

Because the inferred halo flattening is proportional 
to the square of the true width of the gas layer [Eq. 
(||)], the thickness of the gas layer has to be measured 
reliably (Paper II). The difference in the inferred q p 
between the thickness for the fixed-inclination case 
(q p = 0.3to'i) and the free-inclination case (q p = 0.2 
±0.1) exceeds the g p -error resulting from the errors in 
the measured gaseous velocity dispersions by a factor 
of two. 

3.3. The Stellar Mass Distribution 

3.3.1. The vertical distribution 

Since the shape of the DM halo can be determined 
irrespective of the actual mass of the stellar disk, its 
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vertical distribution is unimportant as well. 

3.3.2. Does light trace mass ? 

In Sec. ^] it was assumed that the stellar mass- 
to-light ratio is constant throughout the galaxy. As 
argued above, a vertically varying mass-to-light ratio 
does not affect the inferred halo flattening. 

A radial gradient in the stellar M/C would change 
the radial distribution of stellar mass. Consider the 
simple case where the true radial scale length dif- 
fers from the value used for /ir; this difference would 
expand the horizontal scale of the predicted flaring 
curves by a factor /iR,trueMri, false, while the mea- 
sured flaring remains unchanged. Thus, the usage 
of too long a radial scale length leads to an overes- 
timate of the halo flattening (too round a halo). In 
many early type spirals the scale length as determined 
from near-infrared photometry is significantly shorter 
than the scale length at visual wavelengths (Terndrup 
et al. 1994; Peletier et al. 1994). Extinction can not 
be the sole reason for this radial blueing since such 
color gradients are found in face-on systems as well 
(de Jong & van der Kruit 1994; de Jong 1995). In 
fact, it may be that a radially varying star formation 
history (mean stellar age) in combination with a ra- 
dial metallicity gradient cause the observed color gra- 
dients (Josey & Arimoto 1992; de Jong 1995, Chapter 
4). However, such color gradients are smalj^] for late 
type, low luminosity systems like NGC 4244. Since 
extinction is not very important in NGC 4244 and the 
current star formation is low, the measured blue light 
probably traces mass, so that the inferred DM-halo 
flattening is probably unaffected. 

3.4. Distance Dependence 

By combining the equations (D30) and (B16) from 
Paper I, it follows that the model thickness, inside the 
stellar disk, is proportional to -\//irz , equivalent to a 
linear distance dependence. Likewise, the the thick- 

5 There are twelve systems in de Jong's (1995) sample which 
are comparable to NGC 4244: in scale length, total lu- 
minosity and total H I mass. For those systems I find 
/iK-band/'iB-band = 1-02 ± 0.06, corresponding to small 
color gradients. 



ness of the gas layer beyond the optical disk is lin- 
early proportional to galactocentric radius [Eq. (^)], 
thus depending linearly on distance. As the linear 
scale of the thickness measurements depend linearly 
on distance as well, the flaring gas layer method yields 
distance-independent values for the mass-to-light ra- 
tio of the stellar disk and the shape of the DM halo. 

3.5. The Influence of the Extragalactic 
Radiation Field 

The H I surface density profile of M33 (Corbelli 
et al. 1989) and the major axis profile of NGC 3198 
(van Gorkom 1993) show sudden change in slope at 
column densities below a few times 10 19 cm -2 . Mal- 
oney (1993) models the NGC 3198 case extensively 
and argues that these slope changes are evidence 
for the extragalactic radiation field (EgRF) ionizing 
the neutral hydrogen layer. As as result, the de- 
gree of ionization will increase with height above the 
plane as well, so that width of the neutral gas layer 
(FWHMjj j ) will be smaller than the thickness of the 
column of total hydrogen (FWHM to t), at low column 
densities. Since the low column density regime gener- 
ally coincides with the region where the thickness of 
the H I layer is sensitive to the shape of the DM halo, I 
investigate the effects of the EgRF on the width of the 
gas layer. Inspection of Maloney's Figs. 4b and 7, and 
comparing the densities of the NGC 3198 model with 
the measured densities in NGC 4244, indicates that 
the thickness of NGC 4244's H I layer may be under- 
estimated by ~ 15% at 12 kpc. However, Maloney's 
models were calculated for an EgRF, $j,4, of 10 4 pho- 
tons cm -2 s _1 while the 2-er upper limit is three times 
larger (Vogel et al. 1995). Assuming that the effects 
of the radiation field are similar if the ratio of criti- 
cal density^ (N^ 11 ) to H I surface density is similar, 
I used Maloney's NGC 3198 model to determine the 
dependence of the ratio of observed to total thick- 

In an analytic approximation of the effects of the EgRF, 
Maloney (1993) finds that at total column densities below 
the critical density, the whole hydrogen column is ion- 
ized, while his detailed analysis indicates a rapid change 
of fractional ionization for column densities near N^ 1 
[= 4.45 10 19 yJ§iA FWHMhi cm" 2 , adapted from Mal- 
oney's Eqn (19)]. 



7 



ness (rpwHii) on the ratio of critical to H I surface 
density (rs )Q With this relation it is possible to cor- 
rect the observed thickness for the effects of the ex- 
tragalactic radiation field, for any observed thickness 
and H I column density, and assumed EgRF strength. 
For the radii of interest I find that the thickness of the 
total hydrogen is (3.4,7.4,18.4,23.7,35.4) x $ M % 
larger than the measured thickness of the H I layer 
at 10.5, 11, 11.5, 12, and 12.5 kpc, respectively. For 
small values of the radiation field, $^4 < 2, the thick- 
ness corrections are small so that inferred halo flat- 
tening is not affected. For stronger radiation fields 
(<&i,4= 2-3) we find: q p = 0.5 ± 0.2. In summary, the 
effects of the extragalactic radiation field on the de- 
termination of q p are at most of the same order of 
magnitude as the uncertainties in the gaseous veloc- 
ity dispersion and inclination. 

3.6. Non-Thermal Pressures Gradients 

The analysis presented in this paper assumes the 
vertical gravitational force is balanced by thermal 
pressure gradients only. It is estimated that in the 
solar neighborhood the pressures due to cosmic rays, 
magnetic fields and kinetic gas pressure are of compa- 
rable magnitude (Spitzer 1978, p. 234). The magnetic 
field of the Milky Way has a very long radial scale 
length (Kulkarni & Heiles 1988, p. 146) and could 
provide significant pressure, even beyond the optical 
disk. However, such a stratified magnetic field is un- 
stable (Parker 1966) and may therefore not contribute 
to the vertical pressure balance. Non-thermal radio 
emission indicates that cosmic rays are closely related 
to sites of star formation (Bicay & Helou 1990), and 
do not extend much beyond such regions. With NGC 
4244's radio continuum flux very low (Paper II) it is 
expected that cosmic ray pressure and/or magnetic 
field strength contributes little to the total pressure, 
especially beyond the optical disk. In the Introduc- 
tion I summarized observational evidence against the 



7 With 

rs = N^'/Ehi, and r FWH M = (FWHM H i/FWHM tot ), 
I find that the following relation holds: tfwhm ~ 1 + 
2.57 (l - exp(-(r E /1.35) 2 ' 2 )) ± 0.015 for r E < 3, while 
for larger rs-values rpwHM decreases again, similar to 
Maloney's Fig. 7. 



importance of non-thermal pressures (Rupen 1991). 
At any rate, if non-thermal pressures are important, 
then an even larger DM-halo midplane density is re- 
quired to produce the observed flaring: additional 
non-thermal pressures require an even flatter DM 
halo (see also Pfenniger et al. 1994, and references 
therein) . 

3.7. Combined systematic effects 

Having established that the velocity dispersion el- 
lipsoid of the H I is probably round 



(Sec. 3.1.), I like to stress once more the importance 
of determining the thickness of the gas layer and the 
velocity dispersion accurately. Furthermore, due to 
projection effects, the derived thickness of the gas 
layer is quite sensitive to the exact value of the in- 
clination (for i85° ± 2° ; Paper II). I illustrate this 
by comparing the inferred DM halo flattenings for the 
case that gas layer widths are determined assuming a 
constant inclination of 84<?5, and the case where the 
whole spectral line cube is used to determine thick- 
ness and (varying) inclination simultaneously. If the 
galaxy were exactly edge-on or at an inclination <; 
80° but larger than 60° , the uncertainty in the de- 
rived gas layer width would be much less severe (Pa- 
per II). 

For example, in the varying (fixed) inclination case 
the dispersion has to be decreased by 6.2 (3.2) times 
the errors on the dispersion measurements^] on all 
9 (10) independent points between 8 and 13 kpc to 
be consistent with a round halo (4.1ct and la times 
smaller for a q p =0.5 halo). Similarly, in order to in- 
fer a round halo, the thickness of the gas layer (as 
presented in Fig. [l]) has to change by +3A<jfwhm 
(+2.5ctfwhm) for the varying (fixed) inclination case 
(+2.1<jfwhm and +1(7 fwhm for a q p =0.5 halo). 
Combining these systematic effects, I find that a 
round halo requires a 2.2<r (1-7<t) change for the vary- 
ing (fixed) inclination case and a 1.5<r (0.5ct) change 
for a q p = 0.5 halo. We see that it is highly unlikely 
that the systematic offsets from the measured values 
of <7gas and FWHM Z required for consistency with a 



8 Alternatively, a 1-a systematic dispersion offset can be 
thought of as a velocity dispersion anisotropy of ~ 8%. 
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round halo, or even for a halo of intermediate flatten- 
ing, are chance occurrences. 

Taking the uncertainties due to the extragalactic 



radiation field (Sec. 3.5.) into account, the required 
systematic offsets in dispersion and thickness to in- 
fer a rounder DM halo become smaller. A moder- 
ately flattened halo with q p = 0.5 is consistent with 
the fixed inclination case for all but the smallest $^4- 
values, while a round halo would be inferred if mea- 
sured dispersion values are systematically decreased 
by V 2 times the dispersion errors. However, when us- 
ing the preferred varying inclination case, a system- 
atic change by 0.9a (1.5a) of the measured disper- 
sions and gas layer widths is required in order to infer 
a q p =0.5 (round) DM halo. 

To summarize, assuming only random errors in in- 
clination, velocity dispersion, gas layer width and ex- 
tragalactic radiation field leads to the conclusion that 
the DM halo of NGC 4244 is highly flattened: q p = 0.2 
±0.1. An extragalactic radiation field equal to the 
current 2-a upper limit increases this value to q p = 
0.3 ± 0.2. A velocity dispersion anisotropy of 25% 
(50%) would be required to increase q p to 0.5 (1.0), 
essentially independent of the value of the extragalac- 
tic radiation field. 

4. DISCUSSION 

The gas layer width increases exponentially with 
radius if the potential is dominated by the stellar disk 
(e.g., van der Kruit 1988; Paper I). For this reason, 
most authors fit exponential functions to the observed 
flaring curves and use simple models for the potential 
of the galaxy to derive that the stellar mass-to-light 
ratio is roughly constant with values between 50 and 
100% of the maximum disk value (van der Kruit 1981; 
Rupen 1991; Paper II). A more detailed model of the 
potential, like the one described in Paper I, could be 
used to re-analyse these data sets. However, in order 
to improve upon the current M/L estimates, both 
the vertical stellar distribution, and the gaseous ve- 
locity dispersion need to be accurately determined. 
For all the extra-galactic systems, the flaring measure- 
ments do not extend beyond the optical disks, so that 
no information regarding the shape of the dark halo 
could be obtained. With the edge of the stellar disk of 



the Milky Way at -14 kpc (Robin et al. 1992) some 
of the H I surveys extend into the halo of the Galaxy. 
The interpretation of this data however would be diffi- 
cult because of the uncertainties in the rotation curve 
(e.g. Merrifield 1992) and the gaseous velocity dis- 
persion. 

Dissipationless galaxy formation simulations tend 
to produce triaxial DM halos (e.g., Frenk et al. 1988; 
Warren et al. 1992), with an intrinsic flattening dis- 
tribution peaked at q p = c/a — 0.5 (± 0.15). Adding 
10% baryonic matter (gas) and including gasdynamics 
to these simulations produces more oblate halos while 
q p remains unchanged (Katz & Gunn 1991; Udry & 
Martinet 1994; Dubinski 1994). Observationally, all 
but the flattest shapes have been reported (for a sum- 
mary, see Sec. [l]). Significantly flattened halos sug- 
gest that the process of galaxy formation is more dis- 
sipational, which might occur if a larger fraction of 
the initial density perturbation is baryonic. Since 
the flaring can be determined at lower inclinations 
as well, possibly at inclinations as low as 60° , the 
shape of dark matter halos surrounding many galax- 
ies of different Hubble types can be determined. It 
will then be possible to gauge the significance of the 
highly flattened dark halo of NGC 4244, and place 
more stringent constraints on the process of galaxy 
formation. 

Rich clusters contain ~ (10±5)/i -15 % (by mass) 
hot X-ray emitting gas (Briel et al. 1992; Mushotzky 
et al. 1995, and references therein), with h the nor- 
malized Hubble constant (h = H /(100 km s _1 Mpc -1 )). 
It might thus be useful to explore the process of 
galaxy formation over a larger range of baryon frac- 
tions. With a total mass-to-light ratio of (200— 300)h 
(Bahcall et al. 1995) for these rich clusters, the mass- 
to-light ratio of cluster baryons (M/Cbbx,c\ ) equals 
(32 - 48)/i -0 - 5 . Standard Big Bang nucleosynthe- 
sis (BBN) models (e.g., Walker et al. 1991) limits 
Af/£Bar,BBN to the range (11 — 35)/i _1 . From these 
limits it follows that for h > 0.5 essentially all baryons 
in the universe are located in rich clusters. In the re- 
gion where the cluster and BBN estimates overlap 
(h > 0.2), M/£ B ar ~ (35 ± 6)hT l . The dynami- 
cal mass-to-light ratios of individual galaxies (Broeils 
1992, 1995) range from 4 to 100 h, with 80% of the 
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systems between 6 and 20 h. It is thus possible that 
the dark halos of individual galaxies consist mainly 
of non-baryonic dark matter, but a 100% baryonic 
dark halo is also possible (Gott et al. 1974; Briel 
et al. 1992; Rubin 1993; Bahcall 1995; Sackett 1995). 
It is not clear however where, and in what form, these 
baryons reside in the galactic halos since all plausible 
forms of baryonic dark matter seem to be excluded 
(Hegyi & Olive 1986). 

Pfenniger et al. (1994) reviewed cold, rotationally 
supported, molecular hydrogen as a dark matter can- 
didate. In their model, small high density molecular 
"clumpuscules" form the building blocks of a highly 
clumped interstellar stellar medium (ISM) in which 
star formation is suppressed as a result of frequent 
collisions between the clumpuscules. Their model is 
best developed in the region beyond the optical disk. 
Inside the optical disk the state of the ISM is much 
more complex as a result of increased radiation field 
and metallicity. With the surface brightness of the 
cosmic background radiation equal to ~ 8.5 Lqpc~ 2 , 
the transition between the two regimes is expected 
to occur around R ~ D25/2. The fact that in many 
galaxies the shape of the rotation curve due to the gas 
is similar to the observed rotation curve (Bosma 1981; 
Carignan et al. 1990; Carignan & Puche 1990; Broeils 
1992, Chap. 10; Paper II) could then be explained if 
3 - 10% of the gaseous surface density is in atomic 
form. Such might be expected in the context of clum- 
puscules hypothesis, or other cold-gas dark matter 
models (e.g., Gerhard & Silk 1995). 

Since the self-gravity of the gas layer beyond the 
optical disk strongly affects the flaring (cf. Paper 
I, Sec. 3.3 and Fig. 8), I investigate whether the 
clumpuscules hypothesis is consistent with the flaring 
data presented in this paper. Penny Sackett kindly 
provided the disk-like surface density distribution in- 
ferred from NGC 4244's rotation curve using a Ke- 
plcrian as well as a flat extrapolation beyond the 
last measured point (S to t,Kop and Etot.Fit! see also 
Sackett 1995). The inferred DM surface density is 
affected by the way the rotation curve is extrapo- 
lated: S dm ,kc P ~ 86exp(— R/5.2), and Sdm.FU ~ 
55 exp (— i?/11.4) M pc~ 2 , but is almost indepen- 
dent of the stellar mass-to-light ratio. The dark-to- 



H I surface density rises in exponential fashion (with 
scale length ~ 1.7 ± 0.2 kpc) from ~ 10 at the edge 
of the optical disk to ~ 100 the edge of the H I disk. 
This increase in dark-to-H I ratio is mostly due to the 
declining H I surface densities below column densities 
of - 2 x 10 20 cm" 2 . Carignan & Puche (1990) find 
a similar effect for NGC 7793 where Sdm/^H I m ~ 
creases exponentially for H I column densities below 
5 x 10 20 cm -2 . This contrasts Bosma's (1981) find- 
ing that, for column densities above ~ 10 20 cm -2 , 
the dark-to-H I surface density ratio is approximately 
constant^. 

The thickness the clumpuscule disk can be calcu- 
latedQ when a vertical velocity dispersion is assumed. 
I find that the clumpuscule disk has a thickness equal 
to the H I layer if the velocity dispersion of the dark 
disk is 1.1 (1.3) times larger than crjj j , for SnM.Kep 
(SdMjFu)- With these dispersions the "Keplerian 
dark disk" is close to being stable against radial insta- 
bilities (Q — 0.8 — 1.2 beyond the stellar disk) while 
the "Flat dark disk" is unstable (Q = 0.8 — 0.3) in the 
5-13 kpc range. Here I include a correction factor 
(1 + (FWHM gaSiZ /2.3) « 1.3) due to the thickness 
of the disk (e.g., Pfenniger et al. 1994). The argu- 
ments above suggest that the dark matter may be 
in the form of a flaring, stable, self-gravitating dark 
disk with a velocity dispersion slightly larger than the 
measured dispersion in the H I . It would be very in- 
teresting to obtain flaring curves and dispersion mea- 
surements for other galaxies to further test the disk- 
like dark matter hypothesis. 

If an anisotropic velocity dispersion ellipsoid of the 
H I rather than a flattened DM halo is the correct ex- 
planation for the data presented in this paper, then 
this anisotropy would have implications for the phys- 
ical state of the ISM beyond the optical disk. An 
anisotropy of the velocity dispersion ellipsoid can oc- 

! Note that the dark halo surface densities were calculated 
in three different ways: Sackett inverts an extrapolated 
rotation curve, Bosma uses Nordsieck's (1973) approxima- 
tion to invert the rotation curve, while Carignan & Puche 
calculate the DM surface density from a z-integration 
through their best fit isothermal DM-halo model. 
'The flaring of the dark disk is calculated self-consistently 
using the method described in Paper I (see also Sec. ET). 
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cur if, for example, the collision time between "parti- 
cles" is large. 

In accordance with van der Kruit's (1981) finding 
that the dark halo of NGC 891 cannot be as flat- 
tened as the stellar distribution, the axial ratio of 
NGC 4244's halo (5:1) is smaller than the axial ratio 
of stellar disk (by a factor of ^2 or ~1.4 depending 
on whether the vertically distribution is exponential 
or sech 2 , respectively). It seems that the "conspiracy 
of shapes" (Sackett et al. 1994), i.e. the similarity of 
the shapes of the dark and light distributions, breaks 
down for the flattest systems. Again, more determi- 
nations of the shape of dark halos will make it possible 
to gauge the frequency of highly flattened halos like 
the one of NGC 4244. 

5. CONCLUSIONS 

As in Paper I, the potential of the global mass dis- 
tribution, the self-gravity of the gas included, is used 
to calculate flaring curves for a series of models with 
various stellar mass-to-light ratios and dark halo flat- 
tcnings. Comparing the observed flaring curve with 
the various model curves I find that the flaring mea- 
surements do not constrain the stellar M/C beyond 
the limits imposed by the rotation curve. The shape 
of the halo on the other hand is well determined. 
The dark matter halo of NGC 4244 is significantly 
flattened, q p =0.2t°^ (E9 - E5), albeit less so than 
the stellar disk. The value for q p reported here is 
somewhat smaller than found earlier (Oiling & van 
Gorkom 1995), mainly as a result of a better treat- 
ment of the gaseous self-gravity and the use of the 
measured gaseous velocity dispersion rather than an 
assumed value. The flaring & velocity dispersion mea- 
surements beyond 3 optical scale lengths are consis- 
tent with the dark matter being co-spatial with the 
H I disk, a configuration which is stable against radial 
instabilities. Uncertainty in the determined inclina- 
tion is the largest source of random errors. The as- 
sumption that the vertical velocity dispersion equals 
the measured tangential dispersion is the most impor- 
tant potential systematic error: the measured flaring 
curve is consistent with a round halo if the anisotropy 
of the gaseous velocity dispersion ellipsoid increases 
with radius. In that case the vertical dispersion of 



the gas has to be 50 - 70% of the measured tangential 
velocity dispersion at the last measured point. 
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Galactocentric radius (kpc) 

Fig. 1. — The measured widths (open and filled triangles with error bars, for the fixed and varying 
inclination case described in Paper II, respectively) together with the model widths corresponding to 
different halo flattenings q p are plotted: q p =1.0, 0.7, 0.5, 0.3, 0.2, and 0.1 from top to bottom. The 
maximum-disk (top panel) and the minimum-disk (bottom panel) disk-halo decompositions are presented 
(the structural parameters are given in the caption of Fig. 15 of Paper II). Both the truncation of the 
stellar disk (at Rmax) and the gaseous self-gravity have been included in the model calculation. The 
measured gaseous velocity dispersions (Paper II) were used in the model calculations. The observed 
flaring curve is the average of the north-eastern and south-western sides. Deemed unreliable in the inner 
5 kpc (Paper II), the measurements from this region are not plotted. All plotted points are essentially 
independent measurements. Beyond the optical disk, the thickness of the gas layer is strongly influenced 
by the shape of the dark halo. Using both inclination cases I find: q p = The mass-to-light ratio 

of the stellar disk is not constrained by these flaring measurements. 
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Fig. 2. — The inferred halo flattening depends strongly upon the value of the velocity dispersion. The 
curves in the top panel (the same labeling of the curves applies as in Fig. |l|) are calculated for the "upper 
limit" of the velocity dispersion as found in the literature. Similarly, the lower panel is calculated for the 
"lower limit". 
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Fig. 3. — The radial variation of the average of the velocity dispersion (full line with the error bars) 
as determined in Paper II is presented. The dotted lines are determined from the south-western side, 
the dashed lines from the north-eastern side. For all galactocentric radii, the full line is the average of 
all curves shown. The filled triangles are the velocity dispersions required to make the model widths for 
a round halo model consistent with the flaring measurements (see Fig. |]: i.e., a zz would decrease by 
almost a factor two over 3 kpc). To indicate the importance of simultaneously determining the thickness 
and inclination, I included the required a zz in case an inclination had to be assumed in determining the 
thickness of the gas layer (filled squares). 
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